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Abstract 



Dimerization of a spin-half Heisenberg antiferromagnet on a square lattice is investigated for several possible dimer- 
ized configurations, some of which are shown to have lower ground state energies than the others. In particular, the 
lattice deformations resulting in alternate stronger and weaker couplings along both the principal axes of a square 
lattice are shown to result in a larger gain in magnetic energy. In addition, a 'columnar' configuration is shown to have 
a lower ground state energy and a faster increase in the energy gap parameter than a 'staggered' configuration. The 
inclusion of unexpanded exchange coupling leads to a power law behaviour for the magnetic energy gain and energy 
gap, which is qualitatively different from that reported earlier. Instead of increasing as d x , the two quantities depend 
on 5 as 8 V / |ln<5| . This is true both in the near critical regime (0 < 5 < 0.1) as well as in the far regime (0 < S < 1). It 
is suggested that the unexpanded exchange coupling is as much a source of the logarithmic dependence as a correction 
due to the contribution of umklapp processes. Staggered magnetization is shown to follow the same 5-dependence in 
all the configurations in the small (5-regime, while for < S < 1, it follows the power law 8 X . 



It is known that dimerization lowers the ground state energy of a spin-half isotropic Heisenberg antiferromagnet 
|p]-px|. In other words, the system stands to gain energy by such lattice deformations that render it dimerized with 
alternate weaker and stronger bonds between spins on neighboring sites. On the other hand the lattice distortions 
cost energy and it is the net energy balance that would determine whether the gain in magnetic energy e{8) — e(0) is 
large enough to affect the spin-Pcicrls transition through dimerization. In a phenomenological theory, this is usually 
seen in terms of an exponent showing the dependence of magnetic and elastic energies on the dimerization parameter 
6; where < S < 1. The parameter S describes the extent of lattice deformation, i.e., it gives the displacement of the 
ith atom through Ui — |(— l) l 5. The spin-dimer formation is usually described by the Hamiltonian 



envisaging alternate stronger and weaker exchange bonds J(l + S) and J(l — 5). These bonds can , in fact, be seen 
to result from the ansatz J (a) = — Jll| . Thus when the distance between a pair of spins decreases from a to <z(l — 6), 
the exchange coupling is taken to change from J to approximately J(l + S). 



PACS numbers: 75.10.Jm, 74.65. +n, 75.50.Ee 



I. INTRODUCTION 




(1) 
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Since < S < 1, and since elastic energies go typically as S 2 , therefore if the magnetic energy gain varies with S 
with an exponent less than 2 then in the limit 5 — > 0, the gain would overwhelm the cost. In situations where the 
ground state is amenable to dimerization, the spin-Peierls transition will be unconditional PJT^]. 

Such aspects as these have been studied extensively in Heisenberg antiferromagnetic chains, as summarized in Table 
1. This aspect has also been revealed by experiments on quasi-one dimensional Heisenberg antifcrromagnet CuGe03 




The situation in two-dimensions is a little more involved because of the possibility of frustration due to a competing 
antiferromagnetic second neighbour interaction which can in principle destroy any LRO of the Neel type as well as the 
possibility of dimerization. Much of the study of two-dimensional Heisenberg antiferromagnet has therefore remained 
focused on the destruction of order by frustration Jig] . Moreover, the ground state of a Heisenberg antiferromagnet 
on a square lattice at zero temperature is Neel-ordered and a critical value of spin-lattice coupling is required for the 
gain in magnetic energy to affect a spin-Peierls transition |I^,[l7],[l8| . It is assumed below that the spin-lattice coupling 
is above the threshold, allowing for dimerization of the lattice. The spin configuration is expected to remain Neel-likc 
under dimerization. This is true in the absence of either frustration or quantum fluctuations which lead to a melting 
of the Neel lattice pf. 

The matter of frustration and quantum fluctuation aside, a simple dimerization of a square lattice is interesting in 
its own right because the lattice distortions can take place in more than one way, each one of the possible configura- 
tions giving a different dependence of the ground state energy on the dimerization parameter. 

Figure (1) shows a few such configurations. Fig. (a) describes a columnar configuration caused by one longitudinal 
static (-7T, 0) phonon, in which the nearest neighbour distances along the x-axis are taken to vary alternately as a(l + 5), 
while those along the y-direction remain a. Fig.(b) shows a staggered configuration in which the lattice deformation 
along the x-direction is alternated as in Fig. (a), but the sequence of alternations is itself alternated as one goes along 
the y-direction. It is caused by a (it, ir) phonon with polarization along the x-axis The difference between the 
earlier considerations of this configuration and ours is that we take into account the elongation in the exchange bond 
along the y-direction also, making it dependent upon the dimerization parameter S. While the coupling along the 
x-direction is alternately and , it is uniformly — 7 === along the y-direction. 

In contrast to the configurations (a) and (b), those in Figs.(c), (d) and (e) allow for simultaneous dimerization along 
both x- and y-directions in the plane. The difference between (c) and (d) is the same as that between (a) and (b): 
configuration (c) is columnar and (d) is staggered. The former, caused by two phonons with wavevectors (it, 0) and 
(0,7r), is called plaquette configuration p2| , p0[ . Fig.(e) shows a much studied configuration, caused by a longitudinal 
(7r, 7r) phonon mode. In these three configurations also the exchange couplings in both x— and y— directions are 
(^-dependent. These five configurations of a dimerized square lattice consisting of N spins are therefore characterized 
by the following nearest neighbour interactions. 



Configuration (a) 

Jy (/. 



That is to say, the dimerization is described by the Hamiltonian 

f 



h =jy: 

i,3 
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J'+l 
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Configuration (b) 



J *,* = J(1-A5),A = ±1 



and the Hamiltonian is given by 
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s/N r 

h=jj: 

i,3 



(i + (-if+is) 
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Configuration (c) 



Jx,\ - Jy,\ - (X+XS) — J ( l _ ^ = ±1 



with the Hamiltonian 



h =jy: 



(l + (-l)^) 



(1 + {-1)38) 



Configuration (d) 



T _ J 

Jx > x — (1+A5) 
^S/,A = 



J(l - XS), A = ±1 



v / 5 2 + (l+A«5) 2 

and the Hamiltonian 



J [1 — XS — (1 — ^-)S" 



'N 



Configuration (e) 

J ^ x = J y< x = y/S'Hl+XS)' 

and the Hamiltonian 
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(1 + (-!)'+' 8) 
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P + (1 + (-l)^) 2 



J (l - XS - (1 - ^)<5 2 ) , A = ±l 



' jv 



Si+i,. 



Si.j+i] 



(3) 



(4) 



(5) 



(6) 



ij ^5 2 + (l + (-l)i+i6y 

Some of the exchange couplings in Eqs.(|[^) blow up at 5 — 1. Our analysis will therefore be confined to < 6 < 1. 

We would like to investigate the five configurations in order to see (i) which of them gives the largest gain in 
magnetic energy as the dimerization sets in, and (ii) whether the use of untruncated exchange coupling leads to a 
single power law valid for the entire range of 8. 



A number of methods can be chosen for this purpose. Spin wave theory, either modified through Takahashi con- 
straint of zero magnetization or a Hartree-Fock approximated non-linear theory, is known to give surprisingly good 
results for spin-half Heisenberg antiferromagnet. Or, a spin wave theory in the spinless fermionic representation 
through Jordan- Wigner transformations takes care of fermionic correlations among the s = | spins. Coupled cluster 
method (CCM) has also been extensively, and successfully, used for spin-half Heisenberg antiferromagnet in one and 
two space dimensions. 

The first two methods belong to the class of mean field theories and hence are not expected to be very reliable when 
it comes to determining critical exponents. The coupled cluster method, on the other hand, is a perturbation method 
in which increasingly higher order correlations can, in principle, be incorporated at will, and which has been shown to 
give satisfactory results even in the lower orders of perturbation |2l|j2^ ] . We believe that the coupled cluster method 
must be sufficiently good to see which of the alternative configurations proposed here is favoured once a spin-Peierls 
transition sets in. 
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II. APPLICATION OF THE COUPLED-CLUSTER METHOD 



In the coupled cluster method it is first necessary to define a ket state starting from a model state | (f> >, which in 
our case is the Neel state. The exact ground state | \& > of the system can then be postulated as 

| >= e s | <j) > (7) 

where S is the correlation operator defined for an N particle system as 

s = ]T 5„ (8) 

71 

withs„= ]T s h ,...., in c}cl cj n (9) 



and C'J is the creation operator defined with respect to the model state. The ground state energy can then be found 
as the eigenvalue of the Hamiltonian in the proposed ground state 



He s | <f) >= E g e s \ 4> >. 

Taking inner product with < <\> \ e~ s gives 

E g =<(f>\ e- s He s \ (f> > . 

The product e~~ s He s can be written as a series of nested commutators in the well-known expansion 

e- s He s = H + [H,S] + i [[H, S],S} + (10) 

where in the present case the series terminates after the fourth term. 



It is usually easier to deal with the s = ^ Heisenberg Hamiltonian by applying a rotation of 180° to the up spin 
sublattice; S x — > —S x , S y — > S y and S z — > —S z such that all the spins in the lattice point down. It is also convenient 
to replace the spin operators with Pauli matrices: 5 J = j — x,y,z pl| , p2fl , A general expression for the nearest 
neighbour spin Hamiltonian in 2D is then 

H = -{H ( 2 (^iV + p + vrT+p) + « +P ) , (ii) 

where p is a vector to the four nearest neighbours. Correspondingly, the string operator S n can now be defined as 

S ^ = iE E stx.-u; • • • 4AA (12) 

ii....i„ ji....j„ 

where subscripts i and j distinguish between sites on the two sublattices. We note that for spin half (<J^) 2 = 
(<rj~) 2 = 0. Truncation of the summation up to the desired level gives rise to different schemes of approximation. 
Taking interaction only between the spins on adjacent sites gives the so-called SUB2-2 scheme. Including interactions 
with the second and fourth neighboring sites gives what is termed as SUB2-4 scheme. And taking the previous two 
schemes including interaction among the four adjacent sites give us what has been termed as local SUB4, or LSUB4 
for short. Each one of these approximations accounts for a different order of perturbation calculation, and takes into 
account a different order of inter-particle correlations. It has been noted that LSUB4 is a sufficiently good approxi- 
mation for calculating the ground state properties of a spin-half Heisenberg system [ p2[ . 

Consider a general case: a Hamiltonian which has four different coupling constants for nearest neighbour interactions 
in two space dimensions. It can be written as 

VW/2 

H = — — ^2 ^ [Jx,\ &2i,j ■ &2i+\,j + Jy,\ &i,2j ■ <7i,2j+\] ■ (13) 
i,j A=±l 
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Here i and j are the two components of the site indices on a square lattice. The correlation operators in the LSUB4 
scheme are defined as 



h3 



11 4i- v<j + h 11 < 2j . + „ / n 



i,2j-v 



(14) 



In these equations, the coefficients ai^^etc, are various forms of the coefficient Si x ...i^ ; jv.-.j in the expressions for 
S2n- The ground state energy within the LSUB4 approximation comes out to be 



E a 



-—[J x ,+i (l + 4ai) + J x _i (1 + 46i) + Jj,,+i (1 + ici) + Jy-i (l+4di)] 



(15) 



The coefficients ai, 02, • • •, I are obtained as solutions of a set of coupled nonlinear equations. These equations arise 
from the fact that such matrix elements as < 4> | O e~ He s \ cf> > are all zero when the operator O is any product 
of creation operators, particularly if it is one of the operator products in the correlation operator S above. 



< a i,2j u %,2j+v 



e~ s He s >= 0; 



-s 



< a 2i,j a 2i+v,j a 2i+2v,j a 2i+3v,j e 

< (T i,2j <T i,2j+u (J i,23+2u (J i,j+3v e ~ He >= °; 



He b >= 0: 



v = ±1,±3 
v = ±1,±3 

v = ±l 

v = ±l 



(16) 



where S = S2 + £3 + £4. These equations translate into the following twelve equations for the unknown parameters: 
(o? - 1 - 2a 3 bi - 2f)J X}+1 + (2oi + 2oi6i - 20363) Jx,-i = 
(26! + 20161 - 2a 3 6 3 ) J a ,+i + (6? - 1 - 2 fll 6 3 - 2g) J x ,-i = 
(2a 3 + 2aia 3 ) J 2 , i+ i + (2a 3 - of + 2a 3 6i - f)J x -i= 
(26 3 - 6? + 2ftl 6 3 - g) J Xt+1 + (2b 3 + 26163) J x ,-i = 

(-2a 3 6i + 2ai/ + a 3 b\ - aia 3 6i) J x ,+i + (f - a\ + 26i/ + a 3 .g + 2aia 3 6 3 ) J X) _i = 
[g - 6? + 2ai.g + 63/ + 2a 3 6i6 3 ) J X)+ i + (-2oi6 3 + 26ip + a \b 3 - 016163)^,-1 = 
(cf - 1 - 2c 3 di - 2h)J Vt+ i + (2ci + 2c x di - 2c 3 d 3 ) J y _i = 
(2di + 2cidi - 2c 3 d 3 ) J y<+X + {d\-l- 2 Cl d 3 - 2l)J y _ x = 
(2c 3 + 2cic 3 ) J y>+1 + (2c 3 - cf + 2c 3 di - = 

(2d 3 - d{ + 2 Cl d 3 - l)J v ,+i + {2d 3 + 2d 1 d 3 )J y ,- 1 = 

(— 2csdi + 2c\h + c 3 d\ — c\c 3 dx)J Vt +\ + (h — c\ + 2d\h + c 3 l + 2ac 3 d 3 )J y -i = 
(I - df + 2ciZ + d 3 h + 2c 3 d 1 d 3 )J y , +1 + {-2 Cl d 3 + 2d x l + c\d 3 - cidid 3 ) J v _i = 

Setting all the coupling constants equal reduces the number of equations from twelve to three and yields exactly 
the same equations as obtained by others |2~l]]22[] . The two sets of six equations each independently determines the 
six coefficients contained in each of them. As expected, the equations are symmetric in some coefficients. The twelve 
coefficients are to be evaluated by solving the above coupled equations numerically for each of the configurations 
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separately by substituting appropriate values of J Xi \ and J y .\- 



To be able to calculate the energy gap between the ground and the first excited states, we shall construct the 
excited ket state |^ e > hi term of a linear excitation operator X, which, operating on the ground state l^o >, takes 
the system to an excited state: >= X >= Xe 5 \<j> > . This operator is constructed as a linear combination 
of products of creation operators |2^] 

X = J>„ (17) 

n 

with 

x n = E xh~~L°M,-<- ( 18 ) 
ji — j„ 

The first excited state is obtained by the operator 

x i = Exj°j" ( 19 ) 
j 

where j can be any site of the two sublattices. It is easily seen that the first excitation energy is 

E e = ~ Q + 2ai + 2b x + 2a 3 + 2b 3 ^j (J Xi+1 + J Xt _ x ) . (20) 

The energy gap for a given 8 is A(<5) = E e (5) — \E g {&)\ . We define gap parameter as D(S) — A(S) — A(0). This is the 
energy required to break a dimcrizcd singlet pair for a given 5. 



To be able to calculate a quantity like staggered magnetization the need of defining the bra state arises. In fact 
the bra state is not simply a conjugate of the ket state defined in Eq.(^). The bra ground state wave function < W \ 
corresponding to the ket state 1^ > can be defined as pT| , p2| 

< $ |=< | S 2n e~ s 

where the correlation operator 5 is built wholly out of destruction operators of the Hamiltonian used. The need of 
defining such operators comes from the fact that e~ s is not equal to e s . In our case the correlation operator is 
defines as 

N/2 

S2n = 1 + S2n 

71=1 

with 

52,1 = W E h-ln^h ■ ■ ■ Vin^h ■■■< T 5„> 

ii....i„ ji....j„ 

where i and j indicate vectors in the two sublattices respectively. The first term in S2n ensures orthonormality of the 
bra and ket state; i.e., <^\^> = <4>\^>^<4>\(t>> = \. The bra state coefficients are found by putting the 
matrix elements of the commutator of the Hamiltonian with a string of creation operators in the states ^ and ^ to 
zero. 

< | Se- S [H, a+a+ ■ ■ ■ <J+a+a+ ■ ■ ■ cr+Je 3 ^ >= 0, n = 1, 2, 3, • • ■ (21) 

Equations (|2l]) form a set of coupled linear equations for the bra coefficients S, with the ket state coefficients already 
known. It is to be noted here that the series of nested commutators in e~ s {H, ■ ■ ■ erf ]e s terminates after a finite 
number of terms. 



The correlation operators in the LSUB4 scheme are: 



G 



S2 = [ ai a 2ij a 2i+l,j + b ^ a 2i,j a 2i-l,j + c l a i,2j a i,2j+l + d ^i,2j a i,2j-l 



*,3 



53 = X! [ a 3 a 2i,j a 2i+3,j + b 3 a 2i,j a 2i-3,j + C 3 a i,2j a i,2j+3 + d 3 a i,2j a i,2j-3 



/ n a *+v,i + 3 n ct 2,-, j - + & n ^,2^+. + 1 n ^--^ 



y=0 



i/=0 



i/=0 



i/=0 



(22) 



In these equations, the coefficients ai.&^etc, are various forms of the coefficient S^...^ ; fo....^ in the expressions for 
s 2n - In the LSUB4 scheme, the staggered magnetization, given by 



M z 



where i runs over one sublattice only, becomes 

M z = 1 — a\di — b\b\ — £1303 — 6363 — 2// — 2gg — c\C\ — d\d\ — C3C3 — d 3 d 3 — 2hh — 211 



The bra state coefficients are determined from the following set of simultaneous equations: 

(—3+6161 + a\a\ + 03(13 + 6363 + b\g + b\aif + 2gg + 2//) J Xl +i+ 

(Si + Si 61 - a 3 ai - 2/ai + 2fb\ai + fa 3 b 3 + 3S6361) J x> _i = 

(61 + Mi - 6361 - 2S6i + 2g&iai + g6 3 a 3 + \ja 1 a 3 )J x . + i + 

{-\ +Siai + 6161 + a 3 a 3 + b 3 b 3 + faj+ghai + 2ff + 2gg)J x _ 1 = 

(-Si 61 -6163 + 03 + 0301 - 2/6i + fbj- /6iai + 06163) J Xi+ i+ 
(-Si 63 + a 3 + a 3 oi + /6 3 ai + 2fg) J x _i = 

(-61 03 + 63 + 6 3 ai + g6ia 3 + 2g/)J x ,+i+ 

(-6101 - aia 3 + 63 + 6361 - 2gai + ga^ - g6idi + /aia 3 ) J x _i = 

(-Si + 2/ai + S6 3 )J x ,+i + (-±Si + / + 2/6i) J x ,_i = 

(-±61 + g + 2flroi) J x ,+i + (-61 + 2g6i + /a 3 ) J x ,_i = 

(-5 + didi + C1C1 + S3C3 + d 3 d 3 + d?7 + dxcji + 2ZZ + 2hh)J Vt+ i + 

(Si + cidi — S3C1 — 2hc\ + 2hd\C\ + fec 3 c?3 + \ld 3 d\)J y -\ = 

(di + dici - d 3 di - 2/di + 2/dici + ld 3 c 3 + \hciC 3 )Jy t+ \+ 

(-5 + C1C1 + didi + S3C3 + d 3 d 3 + ftcf + W1C1 + 2ftft + 2U)J y _i = 

(— Cidi — <ii<i3 + S3 + S3C1 — 2/idi + /idi — hd\C\ + ld\d 3 ) J Vi +i + 
(-Sid 3 +S 3 +S 3 di + ftd 3 ci + 2hl)J y -i = 

(-dic 3 + d 3 + d 3 ci + ldic 3 + 2Zft) Jy,+i + 

(— diCi — S1C3 + d 3 + d 3 di — 2/ci + lc{ — Idxd + hc\c 3 )J y -\ = 

(-Si + 2ftCl + Td 3 )Jy t + 1 + (-|Cl + ft + 2hd 1 )Jy- 1 = 

(-irfi +J+ 2l Cl )J Vt+1 + (-di + 2tai + hc 3 )J y -i = 
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III. RESULTS 



The ground state energy and magnetization can now be calculated as functions of the dimcrization parameter 8 . 
Previous calculations have invariably taken spin-spin exchange couplings alternately as J(l±<5), which, as mentioned 
above, can be taken as an expansion of the interaction in Eq.(||) to order 8, implying that the results are valid only in 
the critical regime 8 — > 0. We notice in our calculations that if in Eqs.(g-g|) all the expansions are terminated at the 
order of 8 then the distinction between configurations (a) and (b) disappears. On the other hand, if the expansion 
is taken to one order higher, then there remains no way to distinguish between configurations (c) and (d). We must 
therefore either go to orders beyond 8 2 in the expansion, or retain the interactions in their unexpanded form. We do 
the latter. An added advantage is that the results will then be valid in the limit S — > 1. 



A. Magnetic energy gain 

Our calculations confirm that, like the chain, the ground state energy of all the five configurations decreases with 
8. This is shown in Figures 2, where e(5) — e(0) is plotted against <5 for the proposed configurations. This conclusion 
is not new for some of the configurations in Fig.l p2] , ^8| , ^3| -|2^1 . However, what is significant is that the ground state 
energy goes down with 8 more rapidly for some configurations than others. In fact, Fig. 2 shows that the ^-dependence 
is markedly different for the two types of dimerized configurations: one in which dimcrization takes place only along 
one axis, and the other, in which it occurs along both the directions. The rate of decrease is significantly higher for 
the latter. Also, the columnar configurations lead to a greater gain in magnetic energy than the staggered ones. It 
also shows that the plaquette configuration of Fig. 1(c) is energetically the most favourable state, as also noted earlier 
|^2|p7[ . Particularly in the complete range of 8 (0 < 8 < 1), the plaquette configuration stands out as the most 
preferred one, while there is hardly a discernible difference among the configurations (a), (b) and (d). 

Configuration (e) is peculiar in the sense that 8 = \ is a special point for it; the shorter bond length is symmetric 
about this point, having a minimum value of At this point the distortions give rise to a rectangular lattice with 

sides \/2 and The energy gain increases with 8 up to S = i, and then goes down. 

It is worth pointing out here that the much simpler mean field methods of spin wave theory - either in the bosonic 
representation through Holstein-Primakoff transformations, or in the fermionic representation through Jordan- Wigner 
transformations - yield very similar results. This has been checked by us separately. 

Earlier calculations on the spin-Peierls instability in a 2D system gave varied results on the critical exponents. 
Monte Carlo calculations of Tang and Hirsch |l2| on the Hubbard model in the limit of infinite on-site repulsion U 
found for the cases corresponding to our configurations (a), (b), (c) and (e) that the magnetic energy gain followed 
a simple power law behaviour and increased as 8 2 . Their cases are different from ours in the sense that couplings 
alternated as J(l±5), and were taken constant along the y-direction in the case (b). Feiguin et al Jl8| obtained similar 
results for configurations (a) and (e) in the Schwinger boson representation. Quantum Monte Carlo calculations of 
Katoh and Imada |p6| showed that in chains that are coupled by an antiferromagnetic coupling the exponent of the 
magnetic energy gain in the 8 — > limit is 1. 

Our results are expected to be different from these because instead of J(l ± 5), we take the unapproximated 
exchange coupling J (a) — — . Our CCM calculations show that the gain in magnetic energy does not vary with 8 as a 
simple power law; it varies as j^-jy for all the five configurations in the range < 8 < 0.1 with the exponent v = 1.5. 
In the complete range < 8 < 1 also, they show the same dependence on 8 with v = 1 for the configurations (a) - 
(d). 

It is interesting to note that earlier results show, as summarized in Table 1, that the dimerization of an antiferro- 
magnetic chain also varies as r^-jr, but only in the small 8 regime (the near critical regime). There, the factor of t^j? 

is brought about in the renormalization group calculations as a correction due to umklapp processes |t],^8| . Our CCM 
results, however, show that even in chains this may be the case when the exchange couplings in the dimerized state 
are taken as j^g, instead of the approximated J(l ± 8). We find for chains that the best fit is obtained with yjj-jy in 
the entire range of 8 rather than only in the range of small <5. With the full exchange couplings, the exponent for the 
chain comes out to be v = | for < 8 < 1, and v = 1.3 - 1.6 for < 8 < 0.1. The latter gives a decent comparison 
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with the numbers in Table 1. 



B. The gap parameter 

The S dependence of the energy gap parameter D(S) defined above for the five configurations is shown in Fig. 3, 
showing greater stabilization of the dimerized state with increasing S. We also find that, like the magnetic energy 
gain, the gap parameter D increases with 5 as j^-gj in the small 5 regime for all the five configurations with v = 1.5. 
The configurations (a) - (d) also have the same dependence on S in the entire range of 5 with v = 1. 

The difference between the dimerization of a square lattice along only one direction (Fig. 1(a) and (b)) and along 
both the directions (Fig. 1(c) ,(d) and (e)) is again markedly brought out in Fig. 3. Also the columnar configurations 
again appear as preferred modes of dimerization over the staggered configurations for having higher values of the gap 
parameter in the region of small <5. 



C. Staggered magnetization 

Our CCM calculations in the LSUB4 approximation give staggered magnetization for the un-dimerized square lat- 
tice M(S = 0)= 0.2965, within about 2% of the exact value of 0.303. As dimerization sets in, magnetization decreases 
in all the configurations we have chosen, as shown in Fig. 4, in agreement with the earlier results for configuration 
(a) p5[ ]. This is also the case for the entire range of S (0 < 5 < 1), except in the case of configuration (e) for which 
the magnetization rises again after S — 

The CCM calculations show that for all the five configurations, the magnetization also varies as j^-jy in the small 
S regime with the exponent v = 1.5, exactly as the energy gain and the gap parameter. However in the regime 

< S < 1, M exhibits a simple power law dependence: M ~ 8 X with x between .65 and .75, as shown in Table 2. 
Configuration (e) has a distinctly different behaviour in this regime. 

To summarize, we have studied the spin-Peierls dimerization of a spin-half Heisenberg antiferromagnet on a square 
lattice taking unapproximated exchange couplings based on the ansatz J (a) — — , and assuming that the spin- 
lattice coupling is above the threshold to affect the spin-Peierls transition. We have included different possibilities 
of dimerization. The ground state energy as well as staggered magnetization decrease continuously with increasing 
dimerization for all the proposed configurations. Of the five configurations, those with dimerization taking place 
simultaneously along both the principal square axes have markedly lower ground state energies and magnetization 
than those with dimerization along only one of the axes, in agreement with the result of Lieb and Nachtergaele p7[ . 
Also, those with columnar dimerization have consistently lower energies than those with the staggered dimerization. 
The plaquette configuration stands out as the most favoured mode of dimerization. The energy gap parameter also 
corroborates the above conclusions. It has also been shown that the magnetic energy gain as well as the gap parameter 
and staggered magnetization depend upon the dimerization parameter S as at least in the S — > regime, the 

1 In <5 1 factor coming in without any considerations of umklapp processes being included. 
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Figure captions 



Figure 1: Five configurations for the dimcrization of a square lattice, (a) a columnar configuration caused by a 
longitudinal (ir, 0) static phonon. The nearest neighbour coupling along the horizontal direction alternates between 
J(l — 8) and J(l + (5), while that along the vertical direction remains J. (b) a staggered configuration caused by a 
(tt,tt) static phonon with polarization along x-direction. Like (a), the dimcrization occurs along one direction only, 
but the sequence of alternate couplings itself alternates along the other direction. The coupling along the vertical 
direction is also taken to vary with 8. (c) Dimcrization along both the directions, caused by (tt, 0) and (0, tt) phonons, 
making a plaquette of four nearest neighbour spins, (d) Again dimcrization along both the directions, but taken 
staggered along the vertical direction, (e) Another staggered dimerization that is caused by a longitudinal (tt, tt) 
phonon. Chains are formed with strong bonds. The square lattice deforms to a rectangular lattice with this mode for 
6=1/2. 

Figure 2: The gain in magnetic energy e(8) — e(0) as dimerization sets in with increasing 8 for the five configurations; 
(a) in the range < 8 < 0.1, and (b) in < 8 < 1. 

Figure 3: Dependence of the energy gap parameter D on 8 for the five dimerization configurations; (a) in the range 
< 8 < 0.1, and (b) in < 8 < 1. 

Figure 4: Staggered magnetization varying with i5 for the five dimerization configurations; (a) in the range < 8 < 
0.1, and (b) in < 8 < 1. 
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Table- 1 : Summary of the critical exponents for spin-Peierls transition in a Heisenberg chain determined by various 
methods. 



Method 


Interval 


e(8) - e(0) 


Exponent 




< 8 < 1 


8 X 


x = 4/3 


T? pnOTTTifi li 7i\ 1~i on crrmm ^1 

JL 11U1 J_±±d± JZjCI L1U11 (LlULlU 1— r 


< 8 < 1 


8 X 


x = 1.53 


2-level RG § 


< 8 < 1 
0.05 < <5 < 0.1 
0.4 < 8 < 0.5 


8 X 

8 2 »/\ ln(<5)| 
<5 2 7| ln(<5) 


x = 1.78 
2v = 1.68±g"|| 
2z/ = 1.31 ±0.02 


Excitation spectrum @ 


< 8 < 1 


8 X 


x = 1.36+^ 


Valence bond ID 


8 < 0.05 
8 > 0.05 


<5 2 7|ln( 

J 1 


= 2/3 
cc = 1.36ig;£ 


Finit size scaling 8 


< 8 < 0.1 
< 8 < 1 


8 2v /\ ln(<J)| 
J 11 


i/ = 0.71 ±0.01 
a; = 1.34 ±0.02 


Exact diagonahzation % 


< 8 < 0.1 


<5 2 7 | In(*)| 


i/ = 2/3 


DMRG 


<5 < 0.05 


5" 


.t = 2/3 



Table 2: Exponents obtained by the CCM method for magnetic energy gain, energy gap and magnetization for the 
five dimerized square lattice configurations. The logarithmic power law goes as 8 V / | ln(5) | in the five configurations 
for both energy gain and gap parameter for both small and full i5. While in the staggered magnetization the 
logarithmic law is valid for <5 — > , but it obeys a simple power law in the full dimerization limit; <5 — > 1. Values of v 
are listed in the following table. 



Configuration 


Interval 


e(S) e(0) 


A(<5) - A(0) 


M(0) - M(<5) 


(a) 


< 8 < 0.1 


v = 1.5 


v = 1.5 


f = 1.5 


< 8 < 1 


v = 1.0 


v = 1.0 




(b) 


< 8 < 0.1 


v = 1.5 


v = 1.5 


v = 1.5 


< 8 < 1 


v = 1.0 


z/ = 1.0 


^.Vb 


(c) 


< 8 < 0.1 


v = 1.5 


v = 1.5 


= 1.5 


< 8 < 1 


v = 1.0 


v = 1.0 




(d) 


< 8 < 0.1 


v = 1.5 


= 1.5 


i/ = 1.5 


< 8 < 1 


v = 1.0 


v = 1.0 




(e) 


< 8 < 0.1 


v = 1.5 


^ = 1.5 


Z/ = 1.5 


< 8 < 1 









12 



